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I briefly review some concepts related to coarse-graining methods for the dynamics of soft mat- 
ter systems and argue that such schemes will almost always need to telescope down the physical 
hierarchy of time-scales to a more compressed, but more computationally manageable, separation. 
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The question of how to properly coarse-grain a dynam- 
ical simulation is a very interesting one. I think there is 
no single answer, but want to use an example here from 
our simulations of colloidal hydrodynamics to make some 
points that I believe are of more general relevance. 

Consider a buoyant colloid of mass M c and a radius 
a = \[im in H 2 0. As described in more detail in its 
behaviour is governed by a series of different timescales 
shown in table |H If you are only interested in the be- 
haviour of the colloids, then the two fastest time-scales, 
the solvent collision time t co i and the solvent relaxation 
time Tf, can be ignored as long as they are shorter 
than any other colloidal time-scales. The first physi- 
cally relevant time-scale is the Fokker Planck time-scale 
Tpp w 10 3 over which the colloid loses memory of the 
short-time forces acting on it @ • For the example colloid, 
the next time-scale up is the sonic time t cs » 6.7 x 10 -10 s. 
Then comes the Langevin time tb ~ 2.2 x 10 _7 s that 
measures the exponential decay time of the velocity auto- 
correlation function within the Langevin approximation. 
Interestingly, for colloids this time-scale is artificial and 
does not have direct physical meaning (see appendix of 
[l[). Next up is the kinematic time t v w 10~ 6 s over which 
vorticity diffuses away from the colloid. If your colloid 
moves a significant fraction of its radius within the time 
t v , then the colloid will feel the effects of its own motion 
from a time t v back, and finite Reynolds number (Re) 
effects start to kick in. For that reason, it needs to be 
kept small compared to time-scales of colloidal diffusion 
or advection. The largest time-scale we consider here is 
the diffusion time td ~ 5s. However, if the colloid also 
moves under an external force with a velocity v s , then 
there is an additional time-scale t s = a/v s that measures 
how long it takes to advect over its radius, and you can 
then also define a related Peclet number Pe — t s /rD 
that measures the relative importance of convection over 
diffusion. 

From the Fokker Planck time on up to the diffusion 
time covers 13 orders of magnitude. It is clearly not 
possible to capture all of these in a simulation. Instead, 
what is needed is time-scale separation. As long as the 
time-scales are properly separated, you should still be 
simulating the correct underlying physics. This process 
can visualized in Figure Q] (taken from ref. Q] ) which 
shows an example of how the hierarchy of time-scales is 



10 s- - 

Kf'V- 

10~' 3 s-- \ 

-12 

10 s - - 
10"" 



10 
10"' 



10 

10"%+-, 



- fluid collision 

- fluid de-correlation 



- Fokker-Planck 



10 s- 



Telescoping down 



10 s- - % 



10 s T T B -Browinian (Langevin) 
'0 s+ \ -kinematic 



10 s. 
10"° 



- Fokker-Planck 



-sonic 

\ -kinematic 

t„ -colloid diffusion 



-colloid diffusion 



FIG. 1: Telescoping down: The hierarchy of time-scales for 
a colloid (here the example taken is for a colloid of radius 1 
fim in H2O) is compressed in the coarse-grained simulations 
to a more manageable separation. As long as the physically 
important times are clearly separated, the simulation should 
still generate the correct physical picture. Once the simu- 
lations are completed, they can be related in more detail to 
particular experiments by telescoping back out to the relevant 
experimental time-scales. 



telescoped down to a more computationally manageable 
separation in order maximise simulation efficiency, but in 
such a way that the times are still sufficiently separated 
to correctly resolve the underlying physical behaviour. 
A good example of the rationale behind this thinking 
can be illustrated with the sonic time t cs . Physically it 
needs to be much smaller than the diffusion time, or else 
locally you have supersonic behaviour. But if it is too 
small, the simulation will spend most of its time resolving 
sound waves that may not be that interesting for the 
colloidal behaviour you are trying to reproduce, making 
the simulation very inefficient. 

In order to correctly interpret the physical meaning of 
you simulation you need to telescope the time-scale hi- 
erarchy back out to the physical one you want to study. 
For example, if you are interested in physics that is dom- 
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mated by diffusion, you would map onto your physical 
diffusion time. A consequence of this strategy is that a 
single simulation can map onto many different physical 
times. E.g. if your colloid has a radius a — lpm, then the 
diffusion time is td = 5s, and that fixes the time-scales 
for your simulation. On the other hand, if your colloid 
has a radius a = lOOnm, then td — 5 x 10 _3 s and the ex- 
act same simulation would have a different fundamental 
time-scale. That means that for example the viscosity in 
your simulations would be different depending on what 
colloid size you were mapping to, even though the fun- 
damental physics is the same. This fact suggests that 
physical viscosity is often not such a good parameter to 
try and fit to in a coarse-grained simulation. 

You can also correctly map the same simulation to dif- 
ferent time-scales even though the colloidal system is un- 
changed. Say that you are interested in the longer-time 
behaviour of the velocity autocorrelation function that 
is dominated by the kinematic time t v . In that case if 
you mapped to a physical system with a = 1/im, the 
times would be quite different from what you would get 
when you mapped to the diffusion time of the same sys- 
tem. So here the viscosities etc... would have different 
values depending on what you processes were focusing 
on. This is a good example of a no free lunch theorem. 
If you coarse-grain dynamics, you almost always need to 
do some kind of telescoping down, and that means that 
it is hard to simultaneously match multiple time-scales 
in your system. 

The particular example described here concerns a col- 
loid in suspension where it is relatively straightforward 
to work out what all the time-scales are. Nevertheless, 



we argue in [l[ that many other coarse-graining methods 
for dynamics must make implicit use of the telescoping 
down process. You can make your method work by care- 
fully analyzing your time-scales, and then making sure 
you know how to telescope back out to the experimen- 
tal situation you want to emulate. It helps to do this 
in terms of dimensionless variables. However sometimes 
you can't compress the hierarchy to a computationally 
achievable regime without bringing some time-scales too 
close to each other, or even switching the order of time- 
scales. It is then a matter of subtle judgement if the be- 
haviour that comes out of your simulation is physically 
correct. 

A similar analysis can be used to interpret systems 
where the dynamics are dominated by energy barriers. It 
is going to be very hard to know what the real physical 
time-scales are here because processes depend exponen- 
tially on barrier heights. In a coarse-grained dynamical 
simulation, it may be advantageous to dramatically lower 
(free) energy barriers in order to speed up the simulation. 
The hope is that you keep the relative order that charac- 
terizes the physical system you want to emulate, so that 
at least the qualitative dynamic effects are correct. But 
if, say, you try to interpret your time-scales from a mea- 
surement of single particle diffusion, and then use that 
time to extract a physical time for a different process 
in your simulation that is dominated by energy barri- 
ers, then you will get those numbers completely wrong. 
Clearly, there are many tricky subtleties that arise when 
trying to coarse-grain the dynamics of soft-matter sys- 
tems. 
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TABLE I: Time-scales relevant for colloidal suspension; numerical values are for a buoyant colloid of radius a — lfim in H2O. 
For H2O, the speed of sound is 1.48 x 10 9 fim/s. The kinematic viscosity v = 1 x 10 6 /im 2 /s, and measures the diffusion constant 
with which vorticity diffuses away. 

Solvent time-scales 

Solvent collision time over which solvent molecules interact 
Solvent relaxation time over which solvent velocity correlations decay 
Hydrodynamic time-scales 

Sonic time over which sound propagates one colloidal radius 

Kinematic time over which momentum (vorticity) diffuses one colloidal radius 
Brownian time-scales 

Fokker-Planck time over which force-force correlations decay 

Brownian relaxation time over which colloid velocity correlations decay in Langevin Eq 

Colloid diffusion time over which a colloid diffuses over its radius 
Ordering of time-scales for colloidal particles 

Tcol < Tf,TFP < t cs < TB < T v < TD 
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